Erosion patterns in a sediment layer 
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We report here on a laboratory-scale experiment which reproduces a rich variety of natural pat- 
terns with few control parameters. In particular, we focus on intriguing rhomboid structures often 
found on sandy shores and flats. We show that the standard views based on water surface waves 
do not explain the phenomenon and we evidence a new mechanism based on a mud avalanche 
instability. 
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Many patterns observed 0, Q in natural environments 
stem from erosion/deposition processes. These struc- 
tures are related to transport of solid granular particles 
via a fluid phase that can be either a gas, a liquid, or 
even a flowing granular phase. They span a huge vari- 
ety of spatial and temporal scales. Examples of these 
are fractal river basins [2j, meandering rivers I3J, dune 
fields 3 , granular avalanches Q , ripple marks @ on sand 
banks or on coastal continental platforms. Due in part 
to its economical and environmental impact, elementary 
transport processes involved in erosion are still the focus 
of intense scientific scrutiny. 

It is notoriously difficult to provide a fully consistent 
description of particle laden flows either from a one phase 
or a two phase point of view . Most of the practical 
knowledge on erosion comes from empirical laws often 
derived from field measurements @- Several tentatives 
were made recently to tackle from the statistical physics 
point of view the dynamics of formation of river basins 
[2( but many question are raised when one attempts to 
relate basic transport properties to large scale pattern 
forming instabilities. 

In this letter we report on an experimental setup which 
is designed to produce a generic situation of a falling wa- 
ter level on an erodible sediment layer. This occurs natu- 
rally when the sea retreats from the shore or when a reser- 
voir is drained We use a plexiglass container with 
a flat bottom. A square 130 mmx 130 mm plate of de- 
polished plexiglass slides along the bottom, and is pulled 
by a motor through a translation stage. The whole setup 
can be tilted to an arbitrary angle 9. The first step of the 
experiment is the deposition under water of a sedimented 
powder layer covering the mobile plexiglass sheet. To this 
end, the inclination of the whole setup is lowered so that 
the bottom of the container and the sheet are horizontal 
(8 = 0). A suspension of alumine oxide powder is then 
quickly poured into the container. We use commercial 



abrasive powder made of rough grains with mean diame- 
ter d ~ 30 /xm and density p ~ 2.75 g/cm 3 . The liquid is 
demineralized water, to which we add hydrochloric acid 
so as to lower the pH of the suspension to about 4 in 
order to prevent flocculation of the grains. 
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FIG. 1: Phase diagram localizing the different patterns ob- 
served in the erosion experiments in the tilt angle 9 - velocity 
V space. Lines delimiting the domain boundaries are mere 
guides to the eyes. Letters and symbols correspond to Fig. |2] 
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FIG. 2: Patterns observed in the erosion experiment: a crossed hatched pattern, b disordered branched pattern, c orange skin, 
d chevron structure, e chevrons with oblique channels, f localized pulses at chevron onset. The layer appears darker where it 
has been eroded because the bottom plate is black. A light source to the left creates additional shading. 



height 1 mm, the container and the plate are tilted to 
an angle 9. Part of the liquid is then slowly drained, 
until the contact line of the liquid free surface with the 
tilted bottom of the container reaches about 1cm below 
the top of the mobile sheet. Then the sliding plate is 
pulled out of the liquid at a speed V (see inset of Fig. . 
The layer is filmed by a CCD camera. Its position is 
fixed with respect to the plate and its optical axis is per- 
pendicular to the surface. This simple set-up allows us 
to observe a variety of phenomena and structures de- 
pending on two control parameters, the angle 9 and the 
speed V. The "phase diagram" is sketched on Fig. ^ 
When the sedimented granular layer is pulled out very 
slowly (V < 0.04 cm/s), no pattern is observed. The liq- 
uid seeps out of the sediment, which dries progressively 
without being altered. However, above a critical speed 
and tilt corresponding on Fig. ^ to empty squares, we 
observe the formation of erosive patterns on the granu- 
lar sediment surface. For decreasing angle, the patterns 
faint away progressively and for a tilt angle 9 < 13°, it is 
often difficult to witness of their presence. 

At velocities greater than approx. 0.1 cm/s, surface 
structures appear clearly when the layer is tilted to more 
than about 14°. For angles close to this value, the drain- 
ing liquid leaves a cross-hatched, dense pattern of very 
small and shallow channels (see Fig- EK) ■ 

Around 18°, we obtain a branched, disordered river- 
network, whose biggest branches have widths of about 
1 mm « 30d (Fig-Oa). After the passage of the water con- 
tact line, the surface of the sediment is still smooth and 



the pattern appears with a delay of few tens of seconds. 
The dynamical evolution of the structure can last as long 
as two minutes. First, small localized structures with a 
characteristic angle similar to the previous cross-hatched 
pattern appear almost everywhere and then, bigger and 
bigger disordered structures are created as they merge 
under the action of erosion and sediment transport. We 
will call this regime the disordered regime, refering to the 
random aspect of the final network. Transition between 
cross-hatched and branched is progressive. 

For velocities higher than 0.1 cm/s and slopes increased 
above approximately 19°, there is a sharp (1 — 2°) transi- 
tion to a regime of dimples with a structural aspect simi- 
lar to an orange skin (see Fig-Et)- For steeper slopes, we 
observe the progressive onset of a chevron pattern char- 
acterized by a well defined angle (see Fig.[2i). The cross- 
over region is indicated approximatively by a dashed line 
on the phase diagram. The chevron pattern forms quickly 
(typically five seconds) behind the receding liquid contact 
line. The rhomboid elements characterizing this struc- 
ture have a slightly rounded downhill tip and a height 
profile like fish scales or roof tiles, i.e. the sediment is 
thickest at the downhill tip, with a shallow decrease up- 
hill, and a sharp lower edge. 

Finally, this regime is limited by the maximal stability 
angle 9 m = 35° above which the sediment layer would 
spontaneously avalanche as a whole. 

From the present experiment, we observe no system- 
atic variation of the chevron wavelength A with velocity 
or tilt angle. Nevertheless, the spacing of the chevrons 
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has a tendency to grow as the contact line recedes with 
a mechanism akin to defect fusion; we obtain a mean 
spacing A = 5 mm ± 2 mm. On the other hand, system- 
atic experiments at constant angle 9 = 30°, show clearly 
a decrease of the chevron pattern opening angle ip from 
90° to 30° for increasing velocity (see Fig. |2J. 
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FIG. 3: Chevron alignment angle as a function of velocity. 
Error bars indicate measurement variations 





FIG. 4: Dynamics of the chevron formation, a snapshot of 
the mud avalanche above the contact line, b spatio-temporal 
diagram (split across two lines) showing how chevrons appear 
in the wake of an avalanche. 



Various rhomboid patterns have already been de- 
scribed by geologists in natural environments like on 
sea- ward facing beach slopes in reservoirs or river 
beds after a full drainage of the water. Although there 
seem to be several distinct types, only few attempts at 
explaining a possible mechanism were made Q. Previ- 
ous observations and the few experiments available com- 
monly attribute the formation of chevron patterns to in- 
stabilities at the free surface of the flowing water layer 
(such as hydraulic jumps) which couple to the bottom 
profile [HI ■ Although this could be true for certain types 
of rhomboid ripples in shallow, fast flowing rivers, it can- 
not be the case for the chevron regime described here. 
Indeed, the largest estimation for the Froude number we 
can make at the chevron onset is Fr = V j \fgd = 5 • 10~ 2 , 
which certainly rules out factors like hydraulic jumps of 
the water layer. 

Now we seek to clarify the physical conditions associ- 
ated with the onset of pattern formation. An estima- 
tion of the Darcy flow velocity Vd inside the powder 
yields Vd = gK/v ~ 10~ 5 m/s, which is much smaller 
than the retrieval velocity V. This Vd value is ob- 
tained with a permeability K = 10~ 12 m 2 obtained ex- 
perimentally, and a kinematic viscosity v = 10 _6 m 2 /s. 
Moreover, the capillary length corresponding to a 30 pm 
porous medium under gravity forces is about 20 cm. ft is 
thus legitimate to consider that the sediment remains 
fully soaked with water during retrieval of the plate. 
To estimate the shear exerted by the liquid film at the 
surface, we calculate its thickness h{x,t) in the plate 
reference frame. The retrieval of the plate begins at 
t = 0, and x = is the initial position of the con- 
tact line between liquid surface and sediment. Assum- 
ing that the flow is viscous, so that the average local 
flow velocity is given by V(x,t) — g sin 9 h(x,t) 2 the 



mass conservation equation dth + d x [h[x, t)V(x, £)) = 
yields, for small f ree flow slop es, a self consistent solu- 
tion: h(x,t) = ^J2vx/tg sin 9. Note that we neglected 
both capillary and hydrostatic pressure terms because of 
the small thickness and curvature of the flowing layer. 
This approximation ceases to be valid close to the origin 
and in the vicinity of the junction with the flat water 
level. Note also that the contact line cannot move within 
this approximation, which corresponds to a situation of 
total wetting on the sediment. The maximum height, 
just above the reservoir water level, is thus evaluated to 
be h = \J iy j g sin9) V — 10 /im w d. The first conclu- 
sion is that the Reynolds number, Re = hV/v = 0.03, 
is small enough to justify the lubrication approximation. 
Second, the ratio of the shear exerted by the fluid on 
a grain at the bottom and its apparent weight yields a 
common criterion for the onset of erosion called Shield's 
number, 

= p^tan* /^t) / V\ * = / ApVgtf 

Ap d \V J \ Pw ) v 

p w is the density of water and Ap = p — p w the den- 
sity contrast between grains and liquid. On Fig. [5] we 
plotted the V(6) curve corresponding S — 0.12. The 
scaling implications of this formula should be put to test 
more systematically but so far, it seems to reproduce re- 
markably the shape of the limit were erosion patterns are 
evidenced. Note that a Shields number of value S = 0.12 
is marginally large to represent a situation where a grain 
would be spontaneously dislodged under the action of 
viscous shear. On the other hand, when interpreted 
in the frame work of a Coulomb criterion for the sedi- 
ment layer stability, the shearing strength due to viscous 
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forces could be large enough to trigger an avalanche of 
wet grains. This is indeed what we observe at an an- 
gle above around 18° ± 1°. For this purpose, we visu- 
alized the experiment in the chevron regime by lighting 
the flow through the bottom of the set-up (Fig.QJi): on 
this figure, darker spots corresponding to thicker sedi- 
ment layers. On Fig.^Js we present a space-time view of 
the dynamical process for erosion and chevron formation. 
Each stripe of actual size 0.7 x 6.1 cm 2 corresponds to the 
central part of a picture like 4a. Each successive stripe 
corresponds to a time lag of 0.04s. The line labeled by 
1 indicates the water line retreating at constant veloc- 
ity. The line by labeled 2 shows the front of a sediment 
avalanche triggered at the very top of the soaked layer. 
Clearly this avalanche accelerates as grains accumulate 
at the front, and it eventually catches up with the water- 
line. As the avalanche proceeds, a spatially modulated 
pattern of settling sediment occurs in the wake of the 
surge (line 3 indicates first chevrons) from the top down- 
wards. During this settlement, an horizontal structure 
is created with a well defined spatial selection already 
visible in the vicinity of the avalanche onset. This origi- 
nal structure will give rise to the chevrons pattern since 
the subsequent sedimented structure will occur further 
down with an horizontal spatial shift of roughly one half 
of the wave length. Then a similar scenario is taking 
place for the next structures down and so on. After fi- 
nal settlement of a chevron line, sediment will still flow 
but only between adjacent chevrons. Now, instead of an 
avalanche flow, we rather have a channel flow which dis- 
plays an erosive activity sharpening the pattern until it 
reaches its final form. This final erosion process appears 
as a backwards wave. At larger velocities the channel 
flow is sometimes able to erode the layer down to the 
bottom of the plate, producing a river-like network as in 
Fig.Et. 

There are many open questions left on why such a 
structure is likely to occur in the wake of the avalanche, 
and why it has such a well defined wave-length. Here we 
mention two plausibly relevant paths of thought. Recent 
experiments on sheared dense suspensions in a rotating 
drum show that a particle dense flow is instable under 



shear [llj. Particle segregation occurs in the transverse 
direction, producing regions denser in solid particles that 
may enhance in our situation a localized sedimentary pro- 
cess. Note also that the dynamical aspects related to 
this mud avalanche bears many similarities with studies 
on instabilities and pattern formation in falling viscous 
sheets 12]. To our understanding, this chevron regime 
reveals an interesting and open problem of mud flow and 
sedimentation wave propagation. Along this line, it is 
interesting to note that just at the onset of instability, in 
a tiny velocity region, a regime of localized pulses cross- 
ing each other at 90° is observed (see Fig. |2f), which is 
actually strongly reminiscent of a non-linear waves phe- 
nomenology. 

In conclusion, we presented new results on a labora- 
tory scale experiment producing a rich variety of pat- 
terns due to erosion deposition processes. The notion of 
"phase" we use in the scope of this manuscript is qualita- 
tive (based on visual distinction), and should be defined 
more precisely in terms of an order parameter. We show 
that contrarily to previous propositions these structures 
do not occur as a consequence of water surface waves, 
but we postulate that they are triggered by water vis- 
cous strain in low Reynolds and Froude numbers regime. 
The corresponding Shield's criterion implies certain scal- 
ings of the onset with grain size, viscosity etc., which 
are currently investigated. We especially focussed here 
on an intriguing rhomboid-like (chevron) structure orga- 
nized with a typical distance and a coherent angle be- 
tween chevrons. This avalanche regime is clearly sepa- 
rated from another erosion mechanism producing disor- 
dered networks of branched or cross-hatched channels. In 
the chevron regime we evidence a new mechanism based 
on the triggering of a dense sediment avalanche. The ob- 
served patterns stems from instabilities initiated in the 
wake of the avalanche. Erosion/deposition mechanisms 
further develop in order to carve the pattern to its final 
chevron form. 
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